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^ ■ helicity gravitino production cross sections are calculated by using the full 

. relevant terms in the supergravity lagrangian, and (ii) the high energy photon 
spectrum is obtained by solving the Boltzmann equations numerically. Photo- 

^ dissociation of the light elements (D, T, ^He, ^He) leads to the most stringent 

p • upperbound on the reheating temperature, which is given by (lO^-lO^)GeV 

^ • for the gravitino mass lOOGeV-lTeV. On the other hand, requiring that the 

D . present mass density of photino should be smaller than the critical density, we 

I find that the reheating temperature have to be smaller than (lO^^-lO^'^)GeV 

I for the photino mass (lO-lOO)GeV, irrespectively of the gravitino mass. The 

^ ■ effect of other decay channels is also considered. 
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I. INTRODUCTION 



When one thinks of new physics beyond the standard model, supersymmetry (SUSY) 
is one of the most attractive candidates. Cancellation of quadratic divergences in SUSY 
models naturally explains the stability of the electroweak scale against radiative corrections 
mj^. Furthermore, if we assume the particle contents of the minimal SUSY standard model 
(MSSM), the three gauge coupling constants in the standard model meet at ~ 2 x lO^^GeV 
0,^, which strongly supports the grand unified theory (GUT). 

In spite of these strong motivations, no direct evidence for SUSY (especially superpart- 
ners) has been discovered yet. This means that the SUSY is broken in nature, if it exists. 
Although many efforts have been made to understand the origin of the SUSY breaking, we 
have not understood it yet. Nowadays, many people expect the existence of local SUSY {i.e., 
supergravity) and try to find a mechanism to break it spontaneously in this framework. In 
the broken phase of the supergravity, super-Higgs effect occurs and the gravitino, which 
is the superpartner of graviton, acquires mass by absorbing the Nambu-Goldstone fermion 
associated with the SUSY breaking sector. In a softly broken SUSY models induced by 
minimal supergravity, we expect that the mass of the gravitino ^3/2 lies in the same order 
of those of squarks and sleptons since the following (tree level) super-trace formula among 
the mass matrices J^j's holds [H; 

StrM' ^ E i-l?'i^J + l)trM] ^ 2{n - l)ml/„ (1) 

spin J 

where n is the number of chiral multiplet. For example in models with the minimal kinetic 
term, this is the case and all the SUSY breaking masses of squarks and sleptons are equal 
to gravitino mass at the Planck scale. But contrary to our theoretical interests, we have 
no hope to see the gravitinos directly in the collider experiments since the interaction of 
gravitino is extremely weak. 

On the other hand, if we assume the standard big-bang cosmology, the mass of grav- 
itino is severely constrained. If the gravitino is unstable, it may decay after the big-bang 
nucleosynthesis (BBN) and produces an unacceptable amount of entropy, which conflicts 
with the predictions of BBN. In order to keep success of BBN, the gravitino mass should 
be larger than ~ lOTeV as Weinberg first pointed out 0. Meanwhile, in the case of stable 
gravitino, its mass should be smaller than ~ IkeV not to overdose the universe [0. There- 
fore, the gravitino mass between ~ IkeV and ~ lOTeV conflicts with the standard big-bang 
cosmology. 

However, if the universe went through inflation, we may avoid the above constraints 
P] since the initial abundance of gravitino is diluted by the exponential expansion of the 
universe. But even if the initial gravitinos are diluted, the above problems still potentially 
exist since gravitinos are reproduced by scattering processes off the thermal radiation after 
the universe has been reheated [p|-p!7[]. The number density of secondary gravitino is propor- 
tional to the reheating temperature and hence, upperbound on the reheating temperature 
should be imposed not to overproduce gravitinos. Therefore, even assuming inflation, a de- 
tailed analysis must be done to obtain the upperbound on the reheating temperature. The 
case of stable gravitino has been analyzed in Refs. pB|jr5||T7[ ] and we will not deal with it. 

In this paper, we assume that gravitino is unstable and derive an upperbound on the 
reheating temperature. The analysis of the cosmological regeneration and decay of unsta- 
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ble gravitino has been done in many articles. These previous works show that the most 
stringent upperbound on the reheating temperature comes from the photo-dissociation of 
the hght nuclei (D, T, ^He, ^He). Once gravitinos are produced in the early universe, most 
of them decay after BBN since the lifetime of gravitino with mass O(100GeV — lOTeV) is 
O((10® — 10^) sec). If gravitinos decay radiatively, emitted high energy photons induce cas- 
cade processes and affect the result of BBN. Not to change the abundances of light nuclei, 
we must constrain the number density of the gravitino, and this constraint is translated into 
the upperbound on the reheating temperature. 

In order to analyze the photo-dissociation process, we must calculate the following two 
quantities precisely; the number density of the gravitino produced after the universe re- 
heated, and the high energy photon spectrum induced by radiative decay of gravitino. But 
the previous estimations of these values are incomplete. As for the number density of grav- 



itino, most of the previous works follow the result of Ref . , where the number density is 
underestimated by a factor ~4. Furthermore, in many articles, the spectrum of high energy 
photon, which determines the photo-dissociation rates of light elements, are calculated by 
using a simple fitting formula. In this paper, we treat these effects precisely and found the 
more stringent upperbound on the reheating temperature than the previous calculations. 

The plan of this paper is as follows. In Sec.y, we discuss the interactions of the gravitino. 
In Sec.^, we calculate the gravitino production rate in the early universe. In Sec.[rV|, the 
high energy photon spectrum induced by the radiative decay of gravitino is obtained by 



solving the Boltzmann equations numerically. The results are shown in Sec. VI. Other 



cosmological constraints are considered in Sec. [VIl] and Sec. [VIIl| is devoted to discussions. 



II. INTERACTION OF GRAVITINO 

Before investigating the effects of the gravitino on the inflationary universe, let us discuss 
the interaction of the gravitino briefly. From the supergravity lagrangian we can obtain 
the relevant interaction terms of gravitino iptJ. with gauge multiplets {A^, A) and chiral 
multiplets (0, x)] 

^ = ^^7/. bu, Ip] ^t^F^p + I ■^^^^pL'yu^pXLD^(j)^ + /i.e. I , (2) 

where M = Mpi/y/8n ~ 2.4 x lO^^GeV (with Mpi being the Planck mass).[| Note that other 
interaction terms including gravitino field are not important for our analysis since their 
contributions are suppressed by a factor of M~^. 

Combining Eq.(^ with the renormalizable part of the SUSY lagrangian, we have cal- 
culated the helicity ±| gravitino production cross sections and the results are shown in 
Table |. Note that the cross sections for the processes (B), (F), (G) and (H) are singular 
because of the t-channel exchange of gauge bosons. These singularities should be cut off 
when the effective gauge boson mass Tiieff due to the plasma effect are taken into account. 



^In Ref. 1 11], interactions between the gravitino and the chiral multiplets (the second term in 



Eq.(M)) is ignored in calculating the cross sections for the gravitino production processes. 
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Following Ref. [|TT|, we take 5 = (1 =p cos^)mm = where 5 is a scattering angle in 



the center-of-mass frame, and in our numerical calculations, we choose log(mg^/T^) = 0. 
The effective total cross section in thermal bath Utot is defined by 

Stot = 2 X! ''1^% ^{x+y~^i(>t,+z) , (3) 

where a^^^+y^^ is the cross section for the process x + y ^ ipf^ + Vx = ^ ior incoming 
bosons, r]^ = | for fermions. For the MSSM particle content, T^tot is given by 

^tot = {2.50gliT) + A.mgliT) + ll.W,{T)] , (4) 

where gi, g2 and g^ are the gauge coupling constants of the gauge group U(l)y, SU(2)^ 
and SU(3)^, respectively. Note that in high energy scattering processes, effect of the renor- 
mahzation group flow of the gauge coupling constants should be considered. Using the one 
loop /5-function of MSSM, solution to the renormalization group equation of gauge coupling 
constants is given by 

9^{T)^lg-\M,)-^\og(^)] , (5) 




with bi = 11, 62 = 1, ^3 = ~3. In this paper, we use the above formula. 

From Eq.(^), we can also get the decay rate of the gravitino. In this paper, we only 
consider —>■ J + J decay mode, for which the decay rate is given by 

where m^y is the photino mass. In the case of m;y <^ ^^3/2, this decay rate corresponds to 
the lifetime 

TsMi^, - 7 + 7) = 3.9 X 10« {j^^y sec. (7) 

In the cosmological applications, the gravitino lifetime determines the decay time of the 
gravitino, i.e. the gravitino decays when the Hubble time becomes as the same order of 
Thus the temperature of the background photon at the gravitino decay time becomes lower 
as the lifetime gets longer. 



III. GRAVITINO PRODUCTION IN THE EARLY UNIVERSE 

After the universe has reheated, gravitinos are reproduced by the scattering processes 
of the thermal radiations and decay with the decay rate of order of 771^^2/ -^pi- Since the 
interaction of gravitino is very weak, gravitino cannot be thermalized if the reheating tem- 
perature Tji is less than 0{Mpi). In this case, Boltzmann equation for the gravitino number 
density 71,3/2 can be written as 
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dn3/2 ^ ^rj \ 2 ^3/2 ^3/2 

-—— + 3Hn3/2 = {^totVrel} ^^r■ad " JZ; T ) (8) 

dt (^E3/2) ^3/2 

where H is the Hubble parameter, (■ ■ ■) means thermal average, Urad = C(3)7''^/7r^ represents 
the number density of the scalar boson in thermal bath, Vrei is the relative velocity of the 
scattering radiations {{vrei) = 1 in our case), and ^,3/2/ (-E3/2) is the averaged Lorentz factor. 
Note that the first term of right hand side (r.h.s.) of Eq.(|^) represents contribution from 
the gravitino production process, and the second one comes from the decay of gravitino. 
In Eq.(P), we have omitted the terms which represents the inverse processes since their 
contributions are unimportant at low temperature that we are interested in. In the radiation 
dominated universe, H is given by 



H^^ = J^T^ (9) 
R V 90M2 ' 

where R is the scale factor and A^^, is the total number of effectively massless degrees of 
freedom, respectively. For the particle content of MSSM, A^^,(Tr) ~ 228.75 if Tr is much 
larger than the masses of the superpartners, and N^[T -C IMeV) ~ 3.36. 

At the time right after the end of the reheating of the universe, the first term dominates 
the r.h.s. of Eq.(|D since gravitinos have been diluted by the de Sitter expansion of the 
universe. Using yield variable Y3/2 = n3/2/nrad and ignoring the decay contributions, Eq.(|D 
becomes 

It' ^ Hf ' ^ ^ 

where we have assumed the relation 

RT = const. (11) 
Ignoring the small T-dependence of T^tot, we can solve Eq.([TO|) analytically. Integrating 



Eq.(lO) from the reheating temperature Tr to T (Tr ^ T) and multiplying the dilution 



factor Ns{T)/Ns{Tji), the yield of gravitino is found to be 

For the MSSM particle content, Ns{Tr) ~ 228.75 and NsiT < IMeV) ~ 3.91. Eq.(|g) 
shows that 1^3/2 is proportional to Tr. From Eq.(|l2D, we can derive the simple fitting 
formula for 1^3/2; 

y3,.(r«lMeV). 2.14x10- (^){l-0.02321og(^)}, (13) 
where the logarithmic correction term comes from the renormalization group flow of the 



gauge coupling constants. The difference between the exact formula (|T2D and the above 



approximated one is within ~ 5% (~ 25%) for 10^ GeV <Tr< 10^^ GeV (10^ GeV < Tr < 
10^^ GeV). Note that the numerical value of Y3/2 in our case is about 4-5 times larger than 
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the result in Ref. [|TT| due to the difference of (EtotVrei)- Some comments on this difference 
are given in Sec. VIII . 

As the temperature of the universe drops and approaches the decay term 
becomes the dominant part of the r.h.s. of Eq.(H). Ignoring the scattering term, Eq.(H) can 
be rewritten as 



dt 



3/2 



T3/2 



(14) 



where we have taken 7723/2/ = 1 since gravitinos are almost at rest. Using Eq.(|T^) as 

a boundary condition, we can solve Eq.([T^) and the solution is 



(15) 



riradit) Ns{Tr) 



H{Tj 



R) 



where the relation between t and T can be obtained by solving Eq.(|^) with Eq. (|TT|) ; 



/90M2 



(16) 



IV. RADIATIVE DECAY OF GRAVITINO 



Radiative decay of gravitino may affect BBN. We analyze this effect assuming that the 
gravitino ip^ mainly decays to a photon 7 and a photino 7. 

In order to investigate the photo-dissociation processes, we must know the spectra of the 
high energy photon and electron induced by the gravitino decay. In this section, we will 
derive these spectra by solving the Boltzmann equations numerically. 

Once high energy photons are emitted in the gravitino decay, they induce cascade pro- 
cesses. In order to analyze these processes, we take the following radiative processes into 
account. (I) The high energy photon with energy can become e~ pair by scattering 
off the background photon if the energy of the background photon is larger than m1/e^ 
(with rUe being the electron mass). We call this process double photon pair creation. For 
sufficiently high energy photons, this is the dominant process since the cross section or the 
number density of the target is much larger than other processes. Numerical calculation 
shows that this process determines the shape of the spectrum of the high energy photon 
for Cry > ml/22T. (II) Below the effective threshold of the double photon pair creation, 
high energy photons lose their energy by the photon-photon scattering. But in the limit of 
— 0, the total cross section for the photon-photon scattering is proportional to and this 
process loses its significance. Hence finally, photons are thermalized by (III) pair creation 
in the nuclei, or (IV) Compton scattering off the thermal electron. And (V) emitted high 
energy electrons and positrons lose their energy by the inverse Compton scattering off the 
background photon. Furthermore, (VI) the source of these cascade processes are the high 
energy photons emitted in the decay of gravitinos. Note that we only consider the decay 
channel V'a* — > 7 + 7 and hence the energy of the incoming photon e^o is monochromatic. 

The Boltzmann equations for the photon and electron distribution function and /e 
are given by 
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dt 



dt 



+ 



dt 



+ 



DP 



dt 



+ 



dt 



pp 



cs 



dt 



dt 



PC 



ic 



dt 



+ 



dfe{Ee 



DP 



dt 



+ 



dt 

dfe{Ee 



PC 



9t 



DE 



CS 



9t 



(17) 
(18) 



IC 



where DP, PP, PC, CS, IC, and DE represent double photon pair creation, photon-photon 
scattering, pair creation in nuclei, Compton scattering, inverse Compton scattering, and the 
contribution from the gravitino decay, respectively. Full details are shown in appendix 0. 

In order to see the photon spectrum, we have to solve Eq.(|T^ and Eq.(|TB|). Since the 
decay rate of gravitino is much smaller than the scattering rates of other processes, gravitinos 
can be regarded as a stationary source of high energy photon at each moment. Therefore, 
we only need a stationary solution to Eq.([T7|) and Eq.(0) with non-zero {d / dt)\Y)^ at each 
temperature.0 Note that Eq. fplTf ) and Eq. ([l8|) are linear equations of and /e, and hence, 
once Eq.(p]7D and Eq.(p^) have been solved with some reference value of {df-^/dt)\ji^ we can 
reconstruct the photon spectrum for arbitrary value of {df^/dt)\-£,^ with T and e-yo fixed; 



DE 



{df,/dt)\ 



(19) 



DE 



For each T and e^o, "we have calculated the reference spectra f-y{e^) and fe{Ee) by solving 
Eq.(p!7|) and Eq.(|l8|) numerically with the condition. 



dt 



dfe{Ee 

dt 



0. 



(20) 



Typical spectra are shown in Figs.|I| in which we show the case with e-^o = lOOGeV and lOTeV, 
T = lOOkeV, IkeV, lOeV, and the incoming flux of the high energy photon is normalized to 
be 



df^[e 



t70 



7V^7^ 



dt 



5{e, 



-70, 



GeV^ 



(21) 



DE 



The behaviors of the photon spectra can be understood in the following way. For a given 
temperature T, in the region ^ ml/22T, the photon number density is extremely sup- 
pressed since the rate of double photon pair creation process is very large. Just below this 
threshold value, the shape of the photon spectrum is determined by the photon-photon scat- 
tering process, and if the photon energy is sufficiently small, the Compton scattering with 
the thermal electron is the dominant process for photons. Note that the photon and electron 
spectra are determined almost only by the total amount of energy injection. That is, the 



^This approximation is justified if the scattering rates of high energy photons and electrons are 
sufficiently larger than the expansion rate of the universe. This condition is satisfied in the present 
situation. 
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initial energy dependence of the low energy spectra is negligible. This is consistent with the 
previous work [16|. In Fig.|^, we compare our photon spectrum with the results of the simple 



fitting formula used in Ellis et al. |T6[-Q The fitting formula in Ref. is derived from 
the numerical results given in Refs. 0,|l^ in which, however, the effect of the Compton 
scattering is not taken into account. Our results indicate that the number of Compton scat- 
tering events is comparable to that of the inverse Compton events for such low energy region 
(e-y < ml/ SOT), since the number density of the high energy electron is extremely smaller 
than that of the high energy photon. Therefore, the deformation of the photon spectrum by 
Compton scattering is expected below the threshold of the photon-photon scattering. 



V. BBN AND PHOTO-DISSOCIATION OF LIGHT ELEMENTS 

BBN is one of great successes of the standard big bang cosmology. It is believed that light 
elements of mass number less than 7 are produced at an early stage of the universe when the 
cosmic temperature is between IMeV and lOkeV. Theoretical predictions for abundances of 
light elements are excellently in good agreement with those expected from observations if 
the baryon-to-photon ratio rjB is about 3 x 10^^'^. 

However the presence of gravitino might destroy this success of BBN. Gravitino may 
have three effects on BBN. First the energy density of gravitino at T ~ IMeV speeds up the 
cosmic expansion and leads to increase the n/p ratio and hence ^He abundance. Second, the 
radiative decay of gravitino reduces the baryon-to-photon ratio and results in too baryon- 
poor universe. Third, the high energy photons emitted in the decay of gravitino destroy the 
light elements. Among three effects, photo-dissociation by the high energy photons is the 
most important for gravitino with mass less than ~ ITeV. In the following we consider the 
photo-dissociation of light elements and discuss other effects in Sec. |VIl| . 

The high energy photons emitted in the decay of gravitinos lose their energy during 
multiple electro-magnetic processes described in the previous section. Surviving soft photons 
can destroy the light elements (D, T, ^He, ^He) if their energy are greater than the threshold 
of the photo-dissociation reactions. We consider the photo-dissociation reactions listed in 
Table [I^. For the process D(7,n)p, we use the cross section in analytic form which is given 
in Ref. |^ , and the cross sections for other reactions are taken from the experimental data 
(for references, see Table |I|). We neglect ^He(7, D)D and '^He(7, 2p 2n) since their cross 
sections are small compared with the other reactions. Furthermore, we do not include the 
photo-dissociation processes for ^Li and ^Be because the cross section data for ^Be is not 
available and hence we cannot predict the abundance of ^Li a part of which come from ^Be. 

The time evolution of the light elements are described by 



^Although the photon spectrum is not explicitly given in [16|, we obtain it from their photon 
production spectrum divided by the Compton scattering rate (?i-eO"cs^re/) (where He is the number 
density of the electron and acs the cross section for the Compton scattering process) as Ellis et al. 
did [|^]. Since the cross section has energy dependence, the resultant spectrum (oc e;^'^'^) becomes 
softer than that for photon production (oc e~^'^). 
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^tI] / c?e7^T^a(e7)/7(e7) + I] / f^e^<^T (£7)^^/7 (£7)) (23) 



dt ""^Je, 
dri'Y 

ir~~ -Z^7^^-7"X^aV^7..7V^7. ■ ^ J 
dnSYLe _ 

dt 
dmiie 



i * i * 

n4-a^Y de^^mc-.a{^i)fj{(^j) +Y de.,al^4ue{e^)naft{e^), (25) 

Ul ^ J Ei ^ J Ei 

where cr*^^ is the cross section of the photo-dissociation process i: a-l-7 h+. . . and is the 
threshold energy of reaction i. When the energy of the high energy photon is relatively low, 
i.e. 2MeV ^ < 20MeV the D, T and ^He are destroyed and their abundances decrease. 
On the other hand, if the photons have high energy enough to destroy ^He, it seems that 
such high energy photons only decrease the abundance of all light elements. However since 
D, T and ^He are produced by the photo-dissociation of "^He whose abundance is much 
higher than the other elements, their abundances can increase or decrease depending on the 
number density of high energy photon. When the number density of high energy photons 
with energy greater than ~ 20MeV is extremely high, all light elements are destroyed. But 
as the photon density becomes lower, there is some range of the high energy photon density 
at which the overproduction of D, T and ^He becomes significant. And if the density is 
sufficiently low, the high energy photon does not affect the BBN at all. 

From various observations, the primordial abundances of light elements are estimated 

H 



as 



0.22<y:=(^) < 0.24, (26) 
""^^ >1.8xl0-^ (27) 



p 



^^° + ^^«^^ <1.0xlO-^ (28) 



where pi^ie and ps are the mass densities of ^He and baryon. The abundances of light 
elements modified by gravitino decay must satisfy the observational constraints above. In 
order to make precise predictions for the abundances of light elements, the evolutional 
equations (E^l) - (123) should be incorporated with the nuclear network calculation of BBN. 



Therefore, we have modified Kawano's computer code [|28[ to include the photo-dissociation 
processes. 

From the above arguments it is clear that there are at least three free parameters, i.e. 
mass of gravitino ^3/2, reheating temperature and the baryon-to-photon ratio rjB- Fur- 
thermore we also study the case in which gravitino has other decay channels. In the present 
paper we do not specify other decay channel. Instead, we introduce another free parameter 
-B-y which is the branching ratio for the channel V'/i 7 + 7- Therefore we must study the 
effect of gravitino decay on BBN in four dimensional parameter space. However in the next 
section it will be shown that the baryon-to-photon ratio rjB is not important parameter in 
the present calculation because the allowed value for rjs is almost the same as that in the 
standard case {i.e. without gravitino). 
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VI. RESULTS 



A. Bry = 1 case 

First we investigate the photo-dissociation effect when all gravitinos decay into photons 
and photinos (5^ = 1). We take the range of three free parameters as lOGeV < ^3/2 < 
lOTeV, lO^GeV < Tr < lO^^GeV and 10"^° < r]B < 10"^. In this calculation, we assume 
that the photino is massless. The contours for the critical abundances of the light elements D, 
(D+^He) and ^He in the tib — Tr plane are shown in Figs.0for (a)m3/2 = lOGeV, (b)lOOGeV, 
(c)lTeV and (d)lOTeV, respectively. For low reheating temperature (Tr < lO^GeV), the 
number density of the gravitino is very low and hence the number density of the induced high 
energy photons is too low to affect the BBN. Therefore the resultant abundances of light 
elements are the same as those in the standard BBN. The effect of the photo-dissociation 
due to gravitino decay becomes significant as the reheating temperature increases. 

As seen in Figs.^, the allowed range of the baryon-to-photon ratio is almost the same as 
that without gravitino for 7713/2 ^ ITeV, i.e. very narrow range around t]b ^ 3 x 10"^° is 
allowed. However for 7713/2 ~ ITeV and Tr ~ lO^GeV or 7713/2 ~ ITeV and Tr ~ lO^^GeV, 
lower values of rjB are allowed (Fig.Kc)). In this case, the critical photon energy (~ ml/22T) 
for double photon pair creation process is lower than the threshold of photo-dissociation 
reaction of ^He at the decay time of the gravitino. Therefore, for Tr ^ lO^^GeV, the 
abundance of ^He is not affected by the gravitino decay. Then the emitted photons only 
destroy ^He and D whose abundances would be larger than the observational constraints 
for low baryon density if gravitino did not exist. Therefore one sees the narrow allowed 
band at Tr = lO^GeV where only a small number of '^He and D are destroyed to satisfy 
the constraints (pTl) and (|28|) . For Tr ^ lO^^GeV, since a large number of high energy 
photons are produced even above the threshold of double photon pair creation, a part of 
^He are destroyed to produce ^He and D, which leads to the very narrow allowed region 
at Tr ~ lO^^QgV. However even in this special case, the upper limit of allowed reheating 
temperature changes very little between rjB = 10~^° and 77s ~ 3 x 10^^°. This allows us to 
fix rjB = 3.0 X 10^^° in deriving the upperbound on the reheating temperature. 



The allowed regions that satisfy the observational constraints (|26|)-([28|) also shown in 
Figs.| in the 777,3/2 — Tr plane for r/B = 3 x 10~^°. In Figs.^ and Fig.^(a) one can see four 
typical cases depending on Tr and 7713/2. 

• 7713/2 ^ ITeV, Tr < lO^^GeV: 

In this case the lifetime of the gravitino is so long that the critical energy for the 
double photon process (~ ml/22T) at the decay time of gravitino is higher than the 
threshold of the photo-dissociation reactions for "^He. Thus ^He is destroyed to produce 
T, ^He and D. (Since T becomes ^He by /3-decay, hereafter we mean T and '^He by the 
word "^He".) Since the reheating temperature is not so high, the number density of 
gravitino is not high enough to destroy all the light elements completely. As a result, 
^He and D are produced too much and the abundance of ^He decreases. To avoid 
the overproduction of ^He and D, the reheating temperature should be less than ~ 
(10^ - 109)GeV. 

• 7713/2 < ITeV, Tr > lO^GeV: 

The lifetime is long enough to destroy ^He and the gravitino abundance is very large 
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since the reheating temperature is extremely high. As the result, all the light elements 
are destroyed. This parameter region is strongly excluded by the observation. 

• ITeV < m3/2 ^ 3TeV: 

The lifetime becomes shorter as the mass of gravitino increases, and the decay occurs 
when the double photon pair creation process works well. If the cosmic temperature 
at t = T3/2 is greater than ~ ml/22E4iie (where Ei^^ ~ 20MeV represents the typical 
threshold energy of ^He destruction processes), "^He abundance is almost unaffected 
by the high energy photons as can be seen in Fig.|^(c). In this parameter region, the 
overproduction of (D+^He) cannot occur since ^He is not destroyed. In this case, the 
destruction of D is the most important to set the limit of the reheating temperature. 
This gives the constraint of Tr < (10^ — 10^^)GeV. 

• m3/2 ^ 3TeV: 

In this case the decay occurs so early that all high energy photons are quickly thermal- 
ized by the double photon process before they destroy the light elements. Therefore 
the effect on BBN is negligible. Fig.|^(d) is an example of this case. The resultant 
contours for abundances of light elements are almost identical as those without the 
decay of gravitino. 

B. 5^ < 1 case 

So far we have assumed that all gravitinos decay into photons and photinos. But if other 
superpartners are lighter than the gravitino, the decay channels of gravitino increases and 
the branching ratio for the channel "^/"^ — ^ 7 + 7 becomes less than 1. In this case, vari- 
ous decay products affect the evolution of the universe and BBN. In this paper, instead of 
studying all decay channels, we consider only the 7 + 7 channel with taking the branching 
ratio as another free parameter. With this simplification, the effect of all possible decay 
products other than photon is not taken into account. Therefore the resultant constraints 
on the reheating temperature and the mass of gravitino should be taken as the conserva- 
tive constraints since other decay products may destroy more light elements and make the 
constraints more stringent. 

Although we have four free parameters in the present case, the result for = 1 implies 
that the allowed range of Tr and m3/2 is obtained if we take the baryon-to-photon ratio to 
be 3 X 10^^°. Since our main concern is to set the constraints on Tr and ^3/2, we can safely 
fix 7]B i=3x 10-1°). 

The constraints for B^ = 0.1 and B^ = 0.01 is shown in Fig.^(b) and Fig.^(c) which 
should be compared with Fig.§(a) {B^ = 1 case). Since the number density of the high 
energy photons is proportional to B^, the constraint on the reheating temperature becomes 
less stringent as B^ decreases. In addition, the total lifetime of gravitino is given by 

T-3/2 = T-3/2(V'a. ^ 7 + 7) X ^t- (29) 

Thus the gravitinos decay earlier than that for B^ = 1 case and the constraints from (^He 
+ D) overproduction becomes less stringent. This effect can be seen in Fig.^(b), where the 
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constraint due to the overproduction of (^He + D) has a cut at m^i2 — 400GeV compared 
with ~ ITeV for = 1. 

In Fig.^, the contours for the upperbound on reheating temperature are shown in the 
171^/2 — plane. One can see that the stringent constraint on Tr is imposed for 1713/2 ^ 
lOOGeV even if the branching ratio is small. As mentioned before this constraint should be 
regarded as the conservative one and the actual constraint may become more stringent by 
the effect of other decay products, which will be investigated elsewhere. 



VII. OTHER CONSTRAINTS 

In the previous section, we have considered the constraints from the photo-dissociation 
of light elements. But as we have seen, if the mass of gravitino is larger than a few TeV, 
gravitino decay does not induce light element photo-dissociation and no constraints has 
been obtained. In the case of such a large gravitino mass, we must consider other effects of 
gravitino. 

If we consider the present mass density of photinos produced by the gravitino decay, we 
can get the upperbound on the reheating temperature. In SUSY models with i?-invariance 
(which is the usual assumption), the lightest superparticle (in our case, photino) is sta- 
ble. Thus the photinos produced by the decay of gravitinos survive until today, and they 
contribute to the energy density of the present universe. Since one gravitino produces one 
photino, we can get the present number density of the photino: 

= Y3,2{T < IMeV) x ^T^, (30) 

where Tq is the present temperature of the universe. The density parameter of the 
photino f2;y = m^n^/pc can be easily calculated, where is the photino mass, pc — 
8.1 X IQ-^^/i^GeV^ is the critical density of the universe and h is the Hubble parameter 
in units of lOOkm/sec/Mpc. If we constrain that < 1 in order not to overdose the 
universe, the upperbound on the reheating temperature is given by 

where we have ignored the logarithmic correction term of T^tot- To set the upperbound 
on the reheating temperature, we need to know the mass of photino. If one assumes the 
gaugino-mass unification condition, the lower limit on the mass of photino is 18.4GeV . 
Then we can get the following upperbound on the reheating temperature: 

Tr < 1.5 X lO^^/^^GeV. (32) 

Note that this bound is independent of the gravitino mass and branching ratio. 

Another important constraint comes from the effect on the cosmic expansion at BBN. 
As mentioned before, if the density of gravitino at nucleosynthesis epochs becomes high, the 
expansion of the universe increases, which leads to more abundance of ^He. We study this 
effect by using modified Kawano's code and show the result in Fig.^. In the calculation, we 
take rjB = 2.8 x 10"^° and r„ = 887 sec (where r„ = (889±2.1)sec is the neutron lifetime p4| ) 
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so that the the predicted ^He abundance is minimized without conflicting the observational 
constraints for other hght elements. The resultant upperbound on the reheating temperature 
is given by 

T«<2xl0^3Gev(^)"\ (33) 

for mz/2 > ITeV. 



VIII. DISCUSSIONS 



We have investigated the production and the decay of gravitino, in particular, the effect 
on BBN by high energy photons produced in the decay. We have found that the stringent 
constraints on reheating temperature and mass of gravitino. 

Let us compare our result with those in other literatures. Our number density of gravitino 
produced in the reheating epochs of the inflationary universe is about four times larger than 
that in Ellis et al fll]]. Since Ellis et al. |ll|] note that they have neglected the interaction 



terms between gravitino and chiral multiplets (which is the second term in Eq.(0)), they 
might underestimate the total cross section for the production of gravitino. All previous 
works concerning gravitino problem were based on the gravitino number density given by 
Ellis et al. O]. Therefore our constraints are more stringent than others. In addition, we 



include all standard nuclear reactions as well as photo-dissociation processes in our calcu- 
lation. Therefore the production of (D -|- ^He) contains both contributions from standard 
BBN and photo-dissociation of ^He. Since only the photo-production is taken into account 
in ref. |Tl[], our constraint from (D + ^He) overproduction is more stringent. 

Furthermore, our photon spectrum is different from that in Ref. [0 as shown in Fig.|^.[| 



The spectrum adopted by Ref. [16| has more power to destroy light elements above threshold 
for the photon-photon scattering and less power below the threshold. In Refs. |]18|,|19|, 
Compton scattering process is not taken into account in calculating the photon spectrum 
which Ellis et al. [0 used to derive a fltting formula for the high energy photon spectrum. 
Therefore, it is expected that the difference comes mainly from the neglect of Compton 
scattering off thermal electron, which is the most dominant process for the relatively low 
energy photons (e^ < ml/SOT'^). Our spectrum is also different from that in Ref. [jl4|, 
i.e. our spectrum has larger amplitude especially for heavy gravitino case, and hence our 



constraint on Tr seems to be more stringent .Q The method taken in Ref. is full numerical 
integration (over both time and momentum) of the complicated Boltzmann equations and 
need many time steps to get the flnal spectrum since the typical interaction time is much 
smaller than the cosmic time. Therefore, it may easily contain cumulative numerical errors. 



^However, the upperbound on Tr does not change significantly even if we use the spectrum given 
in Ref. M]- 



^It should be noted that the constraints given in Ref. |14| are obtained by using a simple approx- 
imation and data different from (^ ) - (|28|). 
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Since we obtain the steady state solution at each epoch, there are no cumulative errors in 
the present spectrum. Therefore we believe that the spectrum that we have obtained is more 
precise than those used in other works. 

In summary, we have investigated the photo-dissociation processes of light elements due 
to the high energy photons emitted in the decay of gravitino and set the upperbound on 
the reheating temperature by using precise production rate of gravitino and the spectrum of 
high energy photon. Together with other constraints (the present mass density of photino 
and the enhancement of cosmic expansion due to gravitino) we have obtained the following 
constraint; 

Tr < 10^"^GeV m3/2 < lOOGeV, (34) 
Tr < lO^^^GeV lOOGeV < mg/s < ITeV, (35) 
Tr < lO^-i^GeV ITeV < mg/a ^ 3TeV, (36) 
Tr < lO'^GeV 3TeV < my2 ^ lOTeV. (37) 

This provides a severe constraint in building the inflation models based on supergravity. In 
this paper we have also studied the gravitino which decays into other channels by taking the 
branching ratio as a free parameter. Although this gives conservative upperbound on the 
reheating temperature, the precise constraints cannot be obtained unless various processes 
induced by other decay products are fully taken into account. This will be done in the future 
work pSl . 
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APPENDIX A: BOLTZMANN EQUATION 

In order to calculate the high energy photon spectrum, we must estimate the cascade 
processes induced by the radiative decay of the gravitinos. In our calculation, we have 
taken the following processes into account; (I) double photon pair creation, (II) photon- 
photon scattering, (III) pair creation in nuclei, (IV) Compton scattering off thermal electron, 
(V) inverse Compton scattering off background photon, and (VI) radiative decay of the 
gravitinos. The Boltzmann equations for this cascade processes are given by 



dt 



+ 



dt 



dt 



5/7 



cs 



dt 
dt 



pp 



dt 



PC 



+ 



df- 



IC 



dt 



(Al) 



DE 
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dt 



dt 



+ 



dfe{E, 



DP 



dt 



+ 



PC 



dt 



cs 



dt 



(A2) 



ic 



Below, we see contributions from each processes in detail. 



(I) DOUBLE PHOTON PAIR CREATION [j + -f^e++e-] 

For the high energy photon whose energy is larger than ~ mg/22T, double photon pair 
creation is the most dominant process. 

The total cross section for the double photon pair creation process a dp is given by 



<rnP m = ^rrrl (l - | (s - P') log i±| - 2/3 (2 - p^) 



(A3) 



where Ve = a /me is the classical radius of electron and f3 is the electron (or positron) velocity 
in the center-of-mass frame. Using this formula, one can write down (9/^/9t)|DP as 



dt 



DP 



4m2 



ds sa(j3) 



(A4) 



/3=^l-(4mi/s) 



The spectrum of the final state electron and positron is obtained in Ref. [^], and 
(9/e/9t)|DP is given by 



dfe{Ee 



dt 



-TTr^m^ / ae. 



7 ^3 

7 



de. 







(A5) 



DP " "^■^ v-^ 

where represents the distribution function of the background photon at temperature T, 

7^2 1 



exp(e^/T) — 1 ' 



and function G{Ee, e^, e^) is given by 



(A6) 



log 



4 (^e + K 
EeE'e 
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ml {Ee + ~ ^ 



mlEeE^ 



mz 



Ee + E, 



'\4 



(Ee + i^;;) J i^e'i^;; 



(A7) 



with 



-Eg ~ ^7 + ^7 ~ -^e- 
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(II) PHOTON-PHOTON SCATTERING / 7 + 7^7 + 7 / 
If the photon energy is below the effective threshold of the double photon pair creation, 



photon-photon scattering process becomes significant. This process is analyzed in Ref. [[19 
and for e'^ < 0{ml/T), {df^/dt)\pp is given by 



9 Me' 



dt 



pp 



35584 
101257r^ 

1946 
~506257r 



2^ 2 



2 2 



•J 



(A8) 



For a larger value of e^, we cannot use this formula. But in this energy region, photon-photon 
scattering is not significant because double photon pair creation determines the shape of the 
photon spectrum. Therefore, instead of using the exact formula, we take ml/T as a cutoff 
scale of {df^/dt)\pp, i.e., for < m^/T we use Eg. ( |A8|) and for > ml/T we take 



df,{e' > ml/T) 



dt 



(A9) 



pp 



Note that we have checked the cutoff dependence of spectra is negligible. 



(Ill) PAIR CREATION IN NUCLEI [-/ + N^e+ + e- + N ] 



Scattering off the electric field around nucleon, the high energy photon can produce 
electron positron pair if the photon energy is larger than 2me. Denoting total cross section 
of this process apc, {dfy/dt)\-^p is given by 



dt 



(AlO) 



NP 



where is the nucleon number density. For apc, we use the approximate formula derived 



by Maximon [31 



Differential cross section for this process dapc/dEf, is given in Ref. [^, and {dfe/dt)\]<fp 
is given by 



dt 



UN 



NP 



Ee+m, 



da PC , , 

"^7 J-y[e-yJ 



(All) 



(IV) COMPTON SCATTERING [ j + e' ^ j + e' ] 

Compton scattering is one of the processes by which high energy photons lose their 
energy. Since the photo-dissociation of light elements occurs when the temperature drops 
below ~ O.lMeV, we can consider the thermal electrons to be almost at rest. Using the total 
and differential cross section at the electron rest frame acs and dacs/dEe, one can derive 
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(A12) 
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cs " "'"7 

where fig is the number density of thermal electron. 



(V) INVERSE COMPTON SCATTERING [e^+j^e^+j ] 



Contribution from the inverse Compton process is given by Jones |3^, and {df /dt)\ic is 
given by 
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where function F{e^, E^, e^) is given by 

-^(^7, Eg, e^) = ■! 
with 



2qlogq + {l + 2q){l-q) + ^^{l-q) : for < g < 1, 



: otherwise, 



_ 4e^Ee 
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r,(Ee-e^) 
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(VI) GRAVITINO RADIATIVE DECAY / ^ 7 + 7 / 

Source of the non-thermal photon and electron spectra is radiative decay of gravitino. 
Since gravitinos are almost at rest when they decay and we only consider two body decay 
process, incoming high energy photons have fixed energy e^o; which is given by 



E70 



- m? 



2m 



Therefore, {df^/dt)\-£,^ can be written as 
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FIGURES 



FIG. 1. Typical spectra of photon (the sohd hnes) and electron (the dotted lines). We take 
the temperature of the background photon to be T = lOOkeV, IkeV, lOeV, and the energy of the 
incoming high energy photon e^o is (a) lOOGeV and (b) lOTeV. Normalization of the initial photon 
is given by e^o x {df^{e^)/dt)\i)E = 6{e-y - e^o) GeV^. 

FIG. 2. Photon spectrum derived from the fitting formula used in Ref. [|l^ is compared with 
our result. We take the temperature of the background photon to be lOOeV and the normalization 
of the incoming flux is the same as Fig.^. The solid line is the result of fitting formula, and the 
dotted line is our result with = lOOGeV. 

FIG. 3. Contours for critical abundance of light elements in the r]B — plane for (a) 
m3/2 = lOGeV, (b) mg/g = lOOGeV, (c) mg/a = ITeV and (d)m3/2 = lOTeV. 

FIG. 4. Allowed regions in 771,3/2 — plane for (a) = 1, (h)B^ = 0.1 and (c) = 0.01. In 
the region above the solid curve ^He and D are overproduced, the abundance of ^He is less than 
0.22 above the dotted curve and the abundance of D is less than 1.8 x 10~^ above the dashed curve. 

FIG. 5. Contours for the upper limits of the reheating temperature in the ^3/2 — B^ plane. 
The numbers in the figure denote the limit of the reheating temperature. 

FIG. 6. Upperbound on the reheating temperature. Dashed line represents the constraint from 
the present mass density of photino. Solid curve represents the upperbound requiring ^He < 0.24. 
Constraints from D photo-dissociation is also shown by dotted line. 
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TABLES 



Process 


(A) 




V + A^ 


(B) 


A^ + X^ - 


-^11; + A" 


(C) 


A^^4>i- 


^ + Xi 


(D) 


+ X* - 


V + </'i 


(E) 




V + 


(F) 


X'^X^ - 


-^ip + X" 


(G) 


A'^ + - 


■^i^ + Xj 


(H) 


A'^ + </.i - 


+ j 
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Xi + Xj - 


V + A"^ 


(J) 


+ (pj - 


V + A" 



(5V64vrM2)x 



(8/3) l/"^'^ 

2 



abc 



ji 

rpa 



{-(3/2)+21og(2/5)+5-(l/8)52| 



jabc {-(62/3) + 161og[(2 - 6)/5] + 225 - 26'^ + (2/3)53} 



n~<a 
ji 
2 

rpa 

ji 



(8/3) 
(16/3) 



{-2 + 2\og{2/5)+5] 
{-6 + 81og(2/5) +45- (1/2)52| 

2 

2 



TABLE L Total cross sections for the helicity ib| gravitino production process. Spins of the 
initial states are averaged and those of the final states are summed. /"''^ and Tj" represent the 
structure constants and the generator of the gauge group, respectively. Note that for the processes 
(B), (F), (G) and (H), we cut off the singularities due to the t-, u-channel exchange of gauge bosons, 
taking (1 it 

cos0)mm — 5 where 9 is the scattering angle in the center-of-mass frame. 



Reaction Threshold (MeV) References 



D+7^n + p 2.225 

T+7^n+D 6.257 l23| 

T+7^p + n + n 8.482 

^He +7 p+D 5.494 

^He +7 ^ p+D 7.718 

^He +7 ^ p+ T 19.815 

^He +7 n+^He 20.578 

^He +7 ^ p + n+ D 26.072 



TABLE IL Photo-disociation reactions 
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